Draft version February 1, 2008 

Preprint typeset using LAT^X style emulatcapj v. 04/03/99 



FUSE OBSERVATIONS OF MOLECULAR HYDROGEN IN TRANSLUCENT INTERSTELLAR 
CLOUDS: THE LINE OF SIGHT TOWARD HD 73882 

T. P. Snow 1 , B. L. Rachford 1 , J. Tumlinson 1 , J. M. Shull 1 , D. E. Welty 2 , W. P. Blair 3 , 
R. Ferlet 4 , S. D. Friedman 3 , C. Gry 5,6 , E. B. Jenkins 7 , A. Lecavelier 4 , M. Lemoine 8 , D. 
C. Morton 9 , B. D. Savage 10 , K. R. Sembach 3 , A. Vidal-Madjar 4 , D. G. York 2 , B.-G. 
Andersson 3 , P. D. Feldman 3 , and H. W. Moos 3 

Draft version February 1, 2008 
ABSTRACT 

We report the results of initial FUSE observations of molecular hydrogen (H2) in translucent clouds. 
These clouds have greater optical depth than any of the diffuse clouds previously observed for far-UV H2 
absorption, and provide new insights into the physics and chemistry of such regions. Our initial results 
involve observations of HD 73882, a well-studied southern hemisphere star lying behind substantial 
interstellar material {E B -v = 0.72; Ay = 2.44). We find a total H 2 column density, N(H 2 ) = 1.2 x 10 21 
cm~ 2 , about three times larger than the values for diffuse clouds previously measured in the far-UV. 
The gas kinetic temperature indicated by the ratio N(J=1)/N(J=0) is 58 ± 10 K. With the aid of 
ground-based data to calculate an appropriate multi-component curve of growth, we have determined 
column densities for all rotational levels up to J = 7. The J > 2 states can be reasonably fitted with 
a rotational excitation temperature of 307 ± 23 K. Both the kinetic and rotational temperatures are 
similar to those found in previous investigations of diffuse clouds. The ratios of carbonaceous molecules 
to hydrogen molecules are also similar to ratios in diffuse clouds, suggesting a similar chemistry for this 
line of sight. 

Subject headings: ISM: abundances — ISM: clouds — ISM: lines and bands — ISM: molecules — stars: 
individual (HD 73882) — ultraviolet: ISM 



1. INTRODUCTION 

Molecular hydrogen is the most abundant molecular 
species in the Galactic interstellar medium. While H 2 
is an important constituent of diffuse interstellar clouds 
with visual extinctions Ay less than about 1 magnitude, 
it becomes the dominant form of hydrogen in the so-called 
"translucent" clouds, which are characterized by Ay in 
the range 1-5 mag. H 2 dominates the dynamics, chem- 
istry, and physics of denser clouds, and is the basic raw 
material for star formation. 

Due to the homonuclear structure of H 2 and its lack 
of a dipole moment, ro- vibrational transitions within the 
electronic ground state are quadrupolar with low sponta- 
neous emission coefficients. They arc therefore very diffi- 
cult to observe except where special excitation conditions 
produce detectable ro-vibrational emission in the near-IR 
(e.g. Gautier et al. 1976) or from a space-borne observa- 
tory such as ISO, when a long enough line of sight (Ay « 
20) allows detection of the mid-IR pure rotational emis- 
sion from warm diffuse gas (Falgarone et al. 2000). Cold 
H 2 , the dominant form, can be widely observed through 
electronic transitions in the far UV (A < 1150 A), or with 



difficulty in the IR, where the weak quadrupole transition 
can be observed in absorption only for very high column 
densities (Lacy et al. 1994). 

Previous instruments capable of resolving the far-UV 
absorption bands of H 2 include the Copernicus mission 
(Savage et al. 1977; Spitzer & Jenkins 1975), the IMAPS 
experiment (Jenkins & Peimbert 1997), and the ORFEUS 
spectrometer (Richter et al. 1998; de Boer et al. 1998). All 
of those, however, were limited to relatively bright stars 
with Ay <J 1. 

The FUSE observatory (Moos et al. 2000) is well suited 
for observations of cold H 2 in the diffuse and translucent 
interstellar medium, owing to its high throughput in the 
wavelength region from 905 A to 1187 A, encompassing 
the strong Werner and Lyman bands of H 2 . A survey of 
H 2 in these regions was therefore identified by the FUSE 
PI team as a project of high priority. 

We have assembled a list of 35 lines of sight, which sam- 
ple a wide variety of interstellar environments. These lines 
of sight, with 0.3 mag <; E(B—V) <J 1.1 mag, are also char- 
acterized by a range in extinction properties (determined 
both from Ry and from the shape of the far-UV extinc- 
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tion curve) . Most of the lines of sight have information on 
molecular abundances from previous optical and/or mm- 
wave data. 

In conjunction with the FUSE observing program, we 
and several collaborators are obtaining additional data us- 
ing various ground-based telescopes. We have obtained 
very high-resolution (R ~ 150,000-250,000) spectra of in- 
terstellar K I, Na I, Ca II, and CH absorption along most of 
the lines of sight to understand the cloud velocity structure 
(Welty, Morton, & Snow 2000). We (led by DGY) are ob- 
taining moderately high-resolution (R ~ 50,000), very high 
S/N, nearly complete optical spectra to measure many of 
the diffuse interstellar bands, and the molecules CH, CH + , 
CN, C 2 , and C 3 . T. Oka and B. J. McCall are pursuing 
near-IR measurements of H3 + , and we and F. Chaffee are 
conducting near-IR observations to compare grain mantle 
features such as water ice, the 3.4-/xm hydrocarbon band, 
and the silicate feature at 9.6 /im, with the UV data on 
dust extinction and gas-phase depletions. The ultimate 
goals of this program are to fully understand the masses 
of dense interstellar clouds, to probe the chemistry and 
physics of translucent clouds, and to probe the transition 
region between diffuse and dense interstellar clouds. 

In this Letter, we describe the first observations of H2 
in a translucent cloud line of sight, toward the star HD 
73882. The following sections describe the properties of 
the line of sight deduced from previous observations (§2), 
the FUSE observations and our analysis of the data (§3), 
and the insights gained from these new data (§4). 

2. THE LINE OF SIGHT TOWARD HD 73882 

The star HD 73882 is well known in the translucent 
cloud community because it is of early spectral type 
(08.5V), it is bright enough (V = 7.27) to allow high-S/N 
observations at optical wavelengths, and it has sufficient 
foreground dust (Eb-v = 0.72, Ay = 2.44) and gas to 
have large column densities of many atomic and molecular 
species. Fitzpatrick & Massa (1986, 1988, 1990) derived 
the UV extinction curve and the column density of atomic 
hydrogen from low- resolution IUE spectra. The extinction 
curve shows a steep far-UV rise with significant curvature, 
generally similar to other "dense cloud" curves, as char- 
acterized by Cardelli, Clayton, & Mathis (1988) - a con- 
clusion reached previously for this star by Massa, Savage, 
& Fitzpatrick (1983). However, both the ratio of total to 
selective extinction (Ry — 3.9) and the wavelength of max- 
imum polarization (Serkowski, Mathewson, & Ford 1975) 
suggest that small grains are not as important in this case 
as in most other molecular-cloud lines of sight. The molec- 
ular column densities listed in Table 1 support the view 
that this line of sight is dominated by one or more dense 
clouds, consistent with translucent cloud models (e.g. van 
Dishoeck & Black 1986). 

A scries of papers have provided ground-based data on 
CO, CH, CN, CH+, and C 2 (van Dishoeck et al. 1991; Gre- 
del et al. 1993; see Table 1). The CO emission line pro- 
files reveal three distinct components toward HD 73882, 
at LSR velocities of +5.9, +8.9, and +11.2 km s" 1 . The 
absorption from the other molecular species, however, is 
seen only in a single component near +4.5 km s" 1 . Us- 
ing the 3.6-m reflector and coude echelle spectrometer at 
ESO and the 0.9-m coude feed telescope and coude spec- 
trograph at KPNO, we have obtained very high- resolution 



{R ~ 150,000-250,000) spectra of Na I, K I, and Ca II 
absorption toward HD 73882. A profile-fitting analysis of 
the K I and Na I lines reveals no fewer than 20 velocity 
components, spread over 55 km s _1 . However, the ab- 
sorption is dominated by three narrow (b ~ 0.7 km s _1 ), 
blended components at wlsr ~ +2.7, +4.3, and +5.9 km 
s , with significant contributions from two others at +7.7 
and +11.0 km s _1 . The blend centered at +4.3 km s _1 
probably corresponds to the single (unresolved at R ~ 
60,000 — 100,000) absorption component seen in CH, CN, 
CH + , and C2, and to the CO emission component at +5.9 
km s _1 . The two slightly weaker atomic components may 
correspond to the two other CO components at +8.9 and 
+ 11.2 km s _1 . (Note that the CO column density quoted 
in Table 1 is only for the +5.9 km s _1 component.) The 
numerous weaker features seen in the optical data may 
arise in outlying less dense clouds. 

3. FUSE OBSERVATIONS AND DATA ANALYSIS 

Our FUSE spectrum of HD 73882 is derived from time- 
tagged observations over the course of 8 orbits on 1999 Oct 
30. Several "burst" events occurred during the observation 
(Sahnow et al. 2000). We excluded all photon events that 
occurred during the bursts, reducing effective on-target 
integration time from 16.8 ksec to 16.1 ksec. Strong in- 
terstellar extinction and lack of co-alignment of the SiC 
channels with the LiF channels prevented the collection of 
useful data shortward of 1010 A. 

We performed a simple "collapse" of the 2-dimensional 
spectral image to a 1-dimcnsional spectrum, summing 100 
pixels (LiF 1A segment) or 120 pixels (LiF 2A) in each col- 
umn, and we determined the background level in a similar 
way. This extraction compares favorably to a subsequent 
pipeline processed spectrum. In the regions of interest, 
we obtained a maximum S/N of 20 per 15-pixel resolution 
element in the LiF 1A segment, and 30 in LiF 2A. The 
wavelength solution was derived from narrow lines in the 
diffuse cloud line of sight to LMC star Sk-67 111 (Shull 
et al. 2000). Some residual scatter remained in the wave- 
length solution (Av rms « 6 km s -1 ), precluding a detailed 
analysis of line velocities in this study. Figure 1 depicts a 
portion of our spectrum, along with our model of the H 2 
lines described in the following paragraph. 

The column densities for the J = and J = 1 levels 
come from direct fits to the profiles of the corresponding 
heavily damped lines. The resolution of the FUSE spec- 
tra (R ~ 12,000) is not sufficient to separate the narrow, 
closely spaced components contributing to the interstellar 
line profiles of the high- J lines. In view of the generally 
good correlation between the column densities of Na I and 
H 2 (Federman 1981; Welty & Hobbs 2000) for N(Na I) £ 

10 12 cm' 2 , we have used the five strongest Na I compo- 
nents described in §2 to model the H 2 lines. The velocity 
structure for Na I was preserved for H 2 , and the observed 
b- values for the Na I components were scaled to H 2 by the 
relation 

b= ] /—+vl rb . 

Given the relatively narrow range of column densities, we 
have assumed a one-to-one scaling between Na I and H 2 . 
The predicted H 2 component structure is given in Table 2. 
From this structure we calculated a curve of growth and 
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determined the column densities for J > 2 using 29 H2 
lines. This curve of growth (similar to a single component 
with 6 = 3 km s _1 ) is given in Figure 2, while the column 
densities are given in Table 3. The data in Tables 2 and 3 
were used to calculate the model of the H 2 spectrum given 
in Figure 1. 

Lying on or near the saturation portion of the curve of 
growth, the column densities for J > 4 are highly sensi- 
tive to the assumed component structure. If more than 
five Na I components trace significant amounts of H2, the 
larger "effective" 6-value yields column densities for J = 
4 and J — 5 decreased by an order of magnitude or more. 
Similarly, if the H 2 structure is less complex than we have 
assumed or if components 4 and 5 trace less H 2 (as sug- 
gested by the CH data of Gredel et al. 1993), the opposite 
situation occurs. The J = 5 lines would be shifted ~0.5 
dex to the damping portion of the curve of growth, and 
the column densities for the J = 6 and J = 7 lines could 
be increased by ~0.5 dex or more. 

4. DISCUSSION 

The total H 2 column density of 1.2 x 10 21 cm~ 2 we 
have measured toward HD 73882 is statistically identical 
to the largest previously observed via far-UV absorption 
spectroscopy, 1.1 x 10 21 cm~ 2 toward HD 24534 (Mason, 
et al. 1976). The molecular fraction / = 2N(H 2 )/[2N(H 2 ) 
+ N(H I)] is / = 0.65, nearly identical to that found to- 
ward Zeta Oph (/ = 0.63; Bohlin, et al. 1978). The only 
known larger value, / = 0.80, toward HD 24534 (Snow et 
al. 1998), will be redetermined via a planned FUSE obser- 
vation. 

The abundances of CO, CN, CH, and C 2 , relative to 
the total hydrogen abundance, are more similar to those 
found in diffuse clouds (i.e. toward Zeta Oph) than in 
dark clouds. Because these molecular abundances repre- 
sent only a small fraction of the presumed total carbon 
abundance, we conclude that the clouds toward HD 73882 
have not reached the predicted transition point where car- 
bon becomes primarily molecular (e.g. van Dishoeck & 
Black 1986). The limited data presently available also sug- 
gest that the depletions toward HD 73882 are similar to 
those found in cold diffuse clouds, an issue we will explore 
in more detail in a future paper (Rachford et al. 2000). 

In principle, the CO abundance can be derived from our 
FUSE spectrum. However, the best candidate lines, from 
the C-X (0,0) band at 1088 A, lie on the flat portion of 
the curve of growth where the derived abundance is highly 
sensitive to the chosen 6-value. 

Our derived H 2 column density, combined with the line- 
of-sight extinction properties cited above, show that the 
ratio of hydrogen, both molecular and total, to dust ex- 
tinction, is also similar to the value found for diffuse clouds 
(Bohlin et al. 1978). 

The ratio of molecules in J = 1 to J = (ortho- to 
para-hydrogen) is usually interpreted as a measure of the 
gas kinetic temperature, because in these clouds the col- 
lisional timescale for depopulating these levels is shorter 
than the radiative timescale. Our N(1)/N(0) ratio yields 
Tkin = 58 ± 10 K, consistent with the average value 
2\in = 77 ± 17 K found for diffuse clouds with Coper- 
nicus data (Savage et al. 1977). Shull et al. (2000) find a 
somewhat larger value in a small sample of FUSE data. 

The excitation diagram in Figure 3 shows that the levels 



J > 2 in this line of sight are not in thermal equilibrium at 
Tkin, but instead follow a different distribution character- 
ized by an excitation temperature, T cx = 307 ±23 K. Non- 
thermal excitation of the high- J levels has been explained 
in terms of UV pumping, in which the molecules cascade 
down through upper rotational and vibrational lines fol- 
lowing the absorption of a UV photon (Black & Dclgarno 
1973). In this case, the excited H 2 may arise in the outer, 
optically thin regions of the cloud, rather than in the self- 
shielded interior. If so, the high- J lines may exhibit small 
velocity shifts or enhanced 6-values. However, the limited 
spectral resolution and uncertain wavelength scale of our 
FUSE data precludes a sensitive search for these effects. 

It is noteworthy that even in a line of sight thought to 
be dominated by molecular cloud material within a single 
cloud complex, the velocity structure is very complicated. 
This structure, along with the similarities between the cur- 
rent results and previous results for diffuse clouds, suggests 
the possibility that we looking through several Zeta Oph- 
type clouds rather than a single very dense cloud. 

Our analysis of the excited rotational levels of H 2 de- 
pends heavily on high-resolution ground-based optical ab- 
sorption and mm- wave emission line data to which we had 
access. It is thus imperative for future FUSE interstellar- 
line observations to obtain appropriate high-resolution 
spectra in order to have any hope of unambiguous inter- 
pretation of the lower resolution data. 

We thank the referee, J. H. Black, for helpful comments. 
This work is based on data obtained for the Guaranteed 
Time Team by the NASA-CNES-CSA FUSE mission oper- 
ated by the Johns Hopkins University. Financial support 
to U.S. participants has been provided by NASA contract 
NAS5-32985. 
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Table 1 

Ancillary Abundance Data for HD 73882 



Species 


N 
(cm" 


) 


Reference 


H I 


1.3 


X 


10 21 


Fitzpatrick & Massa 1990 


K I 


2.4 


X 


10 11 


Welty et al. 2000 


Na I 


2.6 


X 


10 13 


Welty et al. 2000 


Ca II 


4.2 


X 


10 12 


Welty et al. 2000 


CH 


3.7 


X 


10 13 


Gredel et al. 1993 


CH+ 


2.4 


X 


10 13 


Gredel et al. 1993 


c 2 


6 


X 


10 13 


Gredel et al. 1993 


CN 


3.8 


X 


10 13 


Gredel et al. 1993 


CO 


3.2 


X 


10 16 


van Dishoeck et al. 1991 



Table 2 

Assumed H 2 component structure 



Comp. 


Rel. strength 


b 


VhSR 






(km s- 1 ) 


(km s- 1 ) 


1 


0.195 


0.98 


+ 2.62 


2 


0.358 


0.98 


+ 4.20 


3 


0.187 


0.98 


+ 5.83 


4 


0.097 


1.21 


+ 7.56 


5 


0.163 


1.16 


+ 10.85 



Table 3 

Column density in each rotational state 



J 


log Nj 

(cm" 2 ) 


Error 

(cm" 2 ) 





20.91 


0.1 


1 


20.59 


0.1 


2 


19.1 


0.2 


3 


18.5 


0.3 


4 


17.5 


0.3 a 


5 


17.0 


0.3 a 


6 


14.5 


0.8 a 


7 


14.4 


0.5 a 


Total 


21.08 


0.1 



a The high- J lines are sub- 
ject to shifts larger than the 
quoted errors with changes in 
the component structure; see 
text. 
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Fig. 1. — Portion of LiF 1A spectrum with H2 model overlaid. The dotted line depicts the continuum derived from modeling the overlapping 
wings of the lines from different bands. The HD line is not included in the model (see Ferlet, et al. 2000 for the HD analysis). 
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Fig. 2. — Multi-component curve of growth for HD 73882. The two curves represent the range in damping constants for the measured lines. 
The error bars lines only give the formal uncertainty in the line fits and may underestimate the total uncertainties for the stronger lines. 
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Fig. 3. — Excitation plot for HD 73882 for J = through J = 7. The solid line corresponds to the assumed kinetic temperature derived 
from J = and J = 1, while the dashed line gives a weighted fit to the J > 2 lines. 



